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ABSTRACT 

It has recently been suggested that the maximum observed quasi-periodic oscillation 
(QPO) frequencies, i^max; for several low-mass X-ray binaries, particularly 4U 1820-30, 
correspond to the orbital frequency at the inner-most stable orbit of the accretion 
disk. This would imply that the neutron stars in these systems have masses ^ 2 Mq, 
considerably larger than any well-measured neutron star mass. We suggest that 
the levelling off of i^qpo may be also understood in terms of a steepening magnetic 
field which, although possibly dipolar at the stellar surface, is altered substantially 
by disk accretion, and presents a "wall" to the accretion flow that may be outside 
the innermost stable orbit. General relativistic effects add to the flattening of the 
'^QPO ~ ^ relation at frequencies below the Kepler frequency at the innermost stable 
orbit. We offer two other possible ways to reconcile the low value of i^max (~ 1060 Hz 
for 4U 1820-30) with a moderate neutron star mass, ~ 1.4M0: at sufficiently large 
M, either (i) the disk terminates in a very thin boundary layer near the neutron star 
surface, or (ii) j^qpo is not the orbital frequency right at the inner edge of the disk, but 
rather at a somewhat larger radius, where the emissivity of the disk peaks. 



Subject headings: accretion, accretion disks - stars: neutron - X-rays: stars 
gravitation - stars: magnetic fields 



Introduction 



Recent observations with Rossi X-ray Timing Explorer (RXTE) have revealed 
kilohertz quasi-periodic oscillations (QPOs) in at least eighteen low-mass X-ray binaries 
(LMXBs; see Van der Klis 1998a, b for a review; also see Eric Ford's QPO web page at 
|http://www. astro. uva.nl/ecford/qp OS. htnii for updated information). These kHz QPOs are 



characterized by their high levels of coherence (with v j up to 100), large rms amplitudes (up 
to 20%), and wide span of frequencies (500 — 1200 Hz). In almost all sources, the X-ray power 
spectra show twin kHz peaks moving up and down in frequency together as a function of photon 
count rate, with the separation frequency roughly constant (The clear exceptions are Sco X-1 and 
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4U 1608-52, van der Klis et al. 1997, Mendez et al. 1998a; see also Psaltis et al. 1998. In Aql 

X-1, only a single QPO has been detected.). Moreover, in several sources, a third, nearly coherent 
QPO has been detected during one or more X-ray bursts, at a frequency approximately equal to 
the frequency difference between the twin peaks or twice that value. (An exception is 4U 1636-53, 
Mendez et al. 1998b.) The observations suggest a generic beat-frequency model where the QPO 
with the higher frequency is associated with the orbital motion at some preferred orbital radius 
around the neutron star, while the lower-frequency QPO results from the beat between the Kepler 
frequency and the neutron star spin frequency. It has been suggested that this preferred radius is 
the magnetosphere radius (Strohmayer et al. 1996) or the sonic radius of the disk accretion flow 
(Miller, Lamb and Psaltis 1998; see also Kluzniak et al. 1990). The recent observational findings 
(e.g., the variable frequency separations for Sco X-1 and 4U 1608-52) indicate that the "beat" is 
not perfect, so perhaps a boundary layer with varying angular frequencies, rather than simply the 
neutron star spin, is involved. 

This paper is motivated by recent RXTE observation of the bright globular cluster source 
4U 1820-30 (Zhang et al. 1998), which has revealed that, as a function of X-ray photon 

count rate, C, the twin QPO frequencies increase roughly linearly for small photon count 
rates {C ~ 1600 — 2500 cps) and become independent of C for larger photon count rates 
(C ~ 2500 — 3200 cps). (The QPOs become unobservable for still higher count rates.) It was 
suggested that the C — independent maximum frequency (fmax = 1060 it 20 Hz) of the upper QPO 
corresponds to the orbital frequency of the disk at the inner-most stable circular orbit (ISO) as 
predicted by general relativity. This would imply that the NS has mass of 2.2Mq (assuming a 
spin frequency of 275 Hz). It has also been noted earlier (Zhang et al. 1997), based on the narrow 
range of the maximal QPO frequencies (t'max ^ 1100 — 1200 Hz) in at least six sources (which 
have very different X-ray luminosities), that these maximum frequencies correspond to the Kepler 
frequency at the ISO, which then implies that the neutron star masses are near 2Mq (see also 
Kaaret et al. 1997). 

The neutron star masses inferred from identifiying fmax with the Kepler frequency at the ISO 
would, if confirmed, be of great importance for constraining the properties of neutron stars and for 
understanding the recycling processes leading to the formation of millisecond pulsars. However, 
while it is tempting to identify t'max with the the orbital frequency at the ISO, this seemingly 
natural interpretation may not be true. One clue that this identification may not be correct is 
that the inferred neutron star masses are substantially above the masses of those neutron stars for 
which accurate determinations are available (Thorsett &; Chakrabarty 1999) even though spin-up 
to Us ~ 300 Hz only requires accretion of a very small amount of material (^ Mq; §2). The 
cause of the flattening of the z^qpo — C correlation, and the value of the maximum frequency, 
are still not understood (and the existence of a plateau in t'Qpo with increasing M is debatable; 
e.g. Mendez et al. 1998c). We suggest in §3 that the steepening of the magnetic field, expected 
near the accreting neutron star, together with general relativistic effect, naturally leads to the 
flattening in the fqpo -M correlation. In §4 we advocate two alternative interpretations of the 
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maximum QPO frequency without invoking excessively large neutron star masses. 



2. Possible Problems with Neutron Star Masses ^ 2 M( 







The most important concern for the inferred neutron star mass of ^ 2Mq is an empirical one. 
LMXBs have long been thought (e.g., Alpar et al. 1982) to be the progenitors of binary millisecond 
radio pulsars. The recent discovery of binary X-ray pulsar SAX J1808-3658 (with spin period 
2.5 ms and orbital period 2 hrs; Wijnands &: van der Klis 1998; Chakrabarty & Morgan 1998) 
appears to confirm this link. Measurements of neutron star masses in radio pulsar binaries give 
values in a narrow range around M 2± IAMq; the data are consistent with a neutron star mass 
function that is flat between ^ I.IM© and ^ I.GMq at 95% CL (Thorsett & Chakrabarty 1999, 
Finn 1994). The masses of neutron stars in X-ray binaries are also consistent with M ~ 1.4 Mq 
(e.g., van Kcrkwijk et al. 1995). Of particular interest is the 5.4 ms recycled pulsar B1855+09 with 
a white dwarf companion: this system is thought to have gone through a LMXB phase (Phinney 
&; Kulkarni 1994), and contains a neutron star with M = 1.41 ibO.lOM© (Thorsett &; Chakrabarty 
1999; earlier Kaspi et al. 1994 estimated M = 1.50 Mq). The 23 ms pulsar PSR B1802-07, 
which is in a white dwarf binary that is also thought to have gone through the LMXB phase, has 
an inferred mass M = 1.26 ±g:^f Mq (95% confidence; Thorsett & Chakrabarty 1999). 

If 1.4 Mq is the mass of the neutron star immediately after its formation in core collapse, 
then to make a 2.2 Mq object would require accretion of material of at least 0.8 Mq. Such large 
accretion mass may be problematic. If we neglect torques on the star due to the interaction of 
its magnetic field and the accretion disk, the added mass needed to spin up the NS to a spin 
frequency Vg = ^s/{2tt) is 



where / = 10^^/45 g cm^ is the moment of inertia, M = I.AMi^^Mq is the neutron star mass and 
and Tin = lOre km is the radius of the inner edge of the accretion disk, which could correspond 
to either the stellar surface (radius R) or the inner-most stable orbit (ISO) in the absence of a 
magnetic field strong enough to influence the flow substantially (see Cook et al. 1994). When the 
neutron star magnetic field is strong enough, the inner radius rin corresponds to the Alfven radius. 
(We note that the positions of all known millisecond pulsars and binary pulsars in the P — P 
diagram for radiopulsars are consistent with spinup via accretion onto neutron stars with dipolar 
surface fields ^ 10^~^ G.) For magnetic accretion, we expect 



where uig = Qs/^Kinn), with O^xirin) the Kepler frequency at rin. The dimensionless function 
f{ujs) includes contributions to the angular momentum transport from magnetic stresses and 
accreting material. It is equal to zero at some equilibrium ujg, but the actual form of /(w^) depends 




(1) 



m, = M^/GM^f{us), 



(2) 
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on details of the magnetic field - disk interaction. Treating constant, we find 
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where, in the last two lines, we have adopted a simple functional form f{tOs) = 1 — ^s', generically. 



AM ^ 0.07 M, 



O" 



(4) 
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where ip{oOs) 1 for ljOs <C Ws,cj assuming that the torque tends to zero at a critical value ojg = ^s,c- 

Large AM is possible if there is a lengthy phase of accretion with nearly zero net torque 
(e.g. accreting O.SMq at a mean accretion rate of M = lO^^gs"^ would require about 400 Myr) 
following a much shorter phase of spin-up to lOg ^s,c (e.g. accreting O.OSMq at M = 10^''gs~^ 
would require 30 Myr). If magnetic field decays during accretion (e.g. Taam & van den Heuvel 
1986, Shibazaki et al. 1989), then the spin-up phase would have been even shorter. (Spin diffusion 
due to alternating or stochastic epsiodes of spin-up and spin-down [e.g. Bildsten et al. 1997, 
Nelson et al. 1997] might be allowed - but constrained - in such a picture.) To accomodate masses 
as large as 2Mq, these LMXBs must be rather old and must have spun up rapidly at first, and 
then not at all for ^ 90% of their lifetimes. Gravitational radiation might provide a mechanism for 
enforcing virtually zero net torque during the bulk of accretion (Bildsten 1998, Andersson et al. 
1999). But equations (|l|) and (^) show that only very small AM is required to achieve Ug ~ 300 Hz, 
irrespective of the mechanism responsible for halting spin-up at such frequencies. 



3. Steepening Magnetic Fields Near the Accreting Neutron Star 

We shall adopt, as a working hypothesis, that the upper QPO frequency is approximately 
equal to the Kepler frequency at a certain critical radius of the disk (Strohmayer et al. 1996; 
Miller, Lamb & Psaltis 1998; van der Klis 1998)Q that is determined by the combined effects of 
general relativity and stellar magnetic field. For sufficiently strong magnetic fields, the disk may 
be truncated near this radius, where matter flows out of the disk and is funneled toward the 
neutron star. This critical radius then corresponds to the usual Alfven radius (Strohmayer et 
al. 1996). Even if the fields are relatively weak (lO'^ — 10^ G) and the field geometry is such that 
matter remains in the disk, the magnetic stress can still slow down the orbital motion in the inner 



^In the model of Titarchuk et al. (1998) , the QPO corresponds of vertical oscillation of the disk boundary layer, 
but the oscillation frequency is equal to the local Kepler frequency. Even in the "non-beat" frequency model of Stella 
and Vietri (1998), the upper QPO frequency still corresponds to the orbital frequency. 
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disk by taking away angular momentum from the flow, and accreting gas then plunges toward 
the star at supersonic speed - a process that is also accelerated by relativistic instability. In this 
case, the critical radius would correspond to the sonic point of the flow (Lai 1998). We neglect the 
possible role of radiative forces discussed by Miller et al. (1998). As emphasized by van der Klis 
(1998a), the fact that similar QPO frequencies (500 — 1200 Hz) are observed in sources with vastly 
different average luminosities (from a few times 10~^LEdd to near LEdd) suggests that radiative 
effects cannot be the only factor that induces the correlation of the QPO frequency and the X-ray 
flux for an individual source. 

Despite many decades of theoretical studies (e.g., Pringle & Rees 1972; Lamb et al. 1973; 
Ghosh & Lamb 1979; Arons 1987; Spruit & Taam 1990; Aly 1991; Sturrock 1991; Shu et 
al. 1994; Lovelace, Romanova &: Bisnovatyi-Kogan 1995, 1999; Miller &: Stone 1997), there 
remain considerable uncertainties on the nature of the stellar magnetic field - disk interactions. 
Among the issues that are understood poorly are the transport of magnetic field in the disk, the 
configuration of the field threading the disk, and the nature of outflows from the disk. To sidestep 
these complicated questions, we adopt a simple phenomenological prescription for the vertical and 
azimuthal components of the magnetic field on the disk. 



where is evaluated at the upper surface of the disk, and (i is the azimuthal pitch angle of the 
field. If we neglect the GR effect, the critical radius r\a is located where the magnetic field stress 
dominates the angular momentum transport in the disk, and it is approximately given by the 
condition 




(5) 




(6) 



using the ansatz equation (|5|), and assuming Keplerian rotation (which may break down in a 
boundary layer near ri^j e.g. LovclciCG ct cil. 1995), wg find. 




(7) 



and the Kepler frequency at rin is 




For a "dipolar" field configuration, n = 3 and i^k {rin) oc as is well-known, but for smaller 

values of n, the dependence steepens; for example, UKirin) oc M for a "monopole" field, n = 2. 
The observed correlation z^qpo oc C may require n < 3, although the relationship between C and 
M is unclear (Mendez et al. 1998c). 



Unusual field topologies are possible as the disk approaches the surface of the neutron star. 
Values of n 7^ 3 (and even violation of power-law scaling) might occur naturally, for open field 
configurations, which may be prevalent because of differential rotation between the star and 
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the disk (e.g. Lovelace et al. 1995). MHD winds driven off a disk could also result in n 7^ 3 
(e.g. Lovelace et al. 1995, Blandford & Payne 1982). Disks that are fully (Aly 1980, Riffert 
1980) or partially (Arons 1993) diamagnetic will also have non-dipolar variation in field strength 
near their inner edges (see also §3.1 below). None of these possibilities requires the field to be 
substantially non-dipolar at the stellar surface, although for disks that penetrate close to the star 
(at Hn — i? ^ R) any non-dipolar field components, if strong enough, would be significant. 

A particular field configuration that could explain the observed variation of i^qpo with C 
might have n < 3 at moderate values of rjn, leading to a strong correlation between z^qpo and M 
(and hence C). As M rises, the disk approaches the star, and the field topology could become 
more complex, resulting in additional, non-power-law radial steepening of the field strength. As is 
argued below, this could happen even if the field is dipolar at the surface of the star, particularly 
if the disk is diamagnetic. This steepening of the field results in a flattening of the i^qpo — M 
relation. Additional flattening results from incipient general relativistic instability at the inner 
edge of the disk. 



3.1. A Specific Ansatz: Diamagnetic Disk 

An illustration of the field steepening discussed above is as follows. Consider a vacuum 
dipole field produced by the star \Bz\ = fi/r"^ (in the equatorial plane perpendicular to the dipole 
axis). Imagine inserting a diamagnetic disk in the equatorial plane with inner radius Flux 
conservation requires TT{rf^ — R'^)\Bz\ = 27rfi/R, which gives the mean vertical field inside between 
R and ri^: 

This field has scaling \Bz\ oc 1/r^ for large r, which would result in UKirm) ^ M (see eq. [^), and 
stiffens as the disk approaches the stellar surface. 

The actual field at r = rjn is difficult to calculate. Aly (1980) found the magnetic field of 
a point dipole in the presence of a thin diamagnetic disk (thickness H <^ r at radius r), and 
demonstrated that the field strength at Hn is enhanced by a factor ~ (rin/-ff)^/^. (See also Riffert 
1980 and Arons 1993.) However, the situation is different for a finite-sized dipole (a conducting 
sphere of radius R) in the presence of a diamagnetic disk. This can be seen by considering 
a simpler problem, where we replace the disk by a diamagnetic sphere (with radius rin)0.The 
magnetic field at radius r (between R and nn) is given by (in spherical coordinates with the 



^ In replacing the disk with a spherical surface, we lose the square- root divergence found by Aly (1980) for 
infinitesmal H/r. But note that for small fields, the disk penetrates near the star, and H may not be very small 
compared with rin — R- In assuming a point dipole, Aly (1980) (and Riffert 1980) exacerbated the divergence, and 
their results probably apply only when rin 3> R- 
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magnetic dipole along the z-axis): 



2/i 2/z \ cos ( 



BAr,0) = + (10) 

where a = R/r\a. Thus the vertical magnetic field at the inner edge of the disk (r = rin) is 

3/i 



BMn)\ = 3 (12) 



m 

We see that the magnetic field steepens as r\a approaches the stellar surface. In reality, some 
magnetic field will penetrate the disk because of turbulence in the disk and Rayleigh- Taylor 
instabilities (Kaisig, Tajima & Lovelace 1992); however, some steepening of the field may remain. 

Adopting the magnetic field ansatz @ and = —^Bz, we can use @ to calculate this 
gives 



2.5 



c 



where Xc = r\^jR and 



/i = BqR^/2 {Bq = W^Bi G is the polar field strength at the neutron star surface), 
Mi.4 = M/{1AMq), Rio = i?/(10km), and Mn = M/{W^'^ gs''^). Alternatively, if we adopt 
dll), we find 

q ^2.5 

Figure 1 shows the Kepler frequency at function the scaled mass accretion rate, 

MuMlf^R'^^^'^Mn/pB^ = O.OT/ft^. Clearly, for smah M, uxinn) depends on M through a 
power-law, but the dependence weakens as M becomes large, in qualitative agreement with the 
observed fqpo-M correlation. General relativistic effects also fiatten the z^qpq-M relation, as we 
discuss next. 



3.2. General Relativistic Effects 



General relativity (GR) introduces two effects on the location of the inner edge of the disk. 
First, the space-time curvature modifies the vacuum dipole field. For example, in Schwarzschild 
metric, the locally measured magnetic field in the equatorial plane is given by 



B' 
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where y = rc^ /{2GM) (Petterson 1974; Wasserman & Shapiro 1983). The GR effect steepens the 
field only at small r. For r = dGM/c^ — lOGM/c^, we find the approximate scaling oc r~^~'', 
with e ~ 0.3 — 0.4. We shall neglect such a small correction to the dipole field given the much 
larger uncertainties associated with the magnetic field - disk interaction. 

A more important effect of GR is that it modifies the the dynamics of the accreting gas 
around the neutron star. Without magnetic field, the inner edge of the disk is given by the 
condition dlx/dr = 0, where Ik is the specific angular momentum of a test mass: 

This would give the usual the ISO at riso = GGM/c^, where no viscosity is necessary to induce 
accretion^. Since magnetic fields take angular momentum out of the disk, we can determine the 
inner edge of the disk using an analogous expression^ 

= -r'B^B,, (18) 

(see eq. j^). In Lai (1998) it was shown that this equation determines the limiting value of the 
sonic point of the accretion flow (although a Newtonian pseudo-potential was used in that paper). 
Adopting the magnetic field ansatz (P), we find 

Similarly, using (12), we have 

It is clear that for 6 1, eq. (|^) or (|20| ) reduces to the Newtonian limit (see §3.1), while for 6 = 
we recover the expected Hn = = QGM/(? . For small 5, the GR effect can modify the inner disk 
radius signficantly. In Fig. 1 we show the orbital frequency at function of the "effective" 

accretion rate) as obtained from (|^) and (pO|). We see that the GR effect induces additional 
flattening in the correlation between VK{^m) and M as r^^ approaches rjso. 



^When viscosity and radial pressure force is taken into account, the flow is transonic, with the sonic point located 
close to Hso. 

*Note that in the limit of perfect conductivity, it is possible to express the Maxwell stress tensor in terms of a 
magnetic field four-vector B that is orthogonal to the fluid velocity four-vector U (e.g. Novikov & Thorne 1973, pp. 
366-367). The field components and Bz in eq. ( |l8| ) and below are actually the projections of B onto a local 
orthonormal basis (i.e. B^ ^ ■ B and Bz ^ Sz ■ B) even though we have retained the nonrelativistic notation for 
these field components. No additional relativistic corrections are required with these identifications understood. 
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We emphasize the phenomenological nature of eqs. (p!^)-(pO|): they are not derived from a 
self-consistent MHD calculation, and take account of the dynamics of the disk under a prescribed 
magnetic field configuration. However, we believe that they indicate the combined effects of 
dynamically altered magnetic field and GR on the inner region of the accretion disk. By measuring 
the correlation between the QPO frequency and the mass accretion rate, one might be able 
to constrain the magnetic field structure in accreting neutron stars, and reach quantitative 
conclusions about the nature of the interaction of the accretion disk and magnetic field. 



4. Where are the QPOs Produced? 

Implicit in the discussion of magnetic fields and fqpo in the preceding sections were the 
assumptions that the QPO arises at a radius outside the star that coincides with the inner radius 
of the accretion disk. Here, we examine two ways in which these assumptions might be violated, 
and show how the relatively small measured values of fmax might be consistent with neutron star 
masses near IAMq. 



4.1. Disk Termination at the Neutron Star Surface 

For the model discussed in §3, the steepening magnetic field and general relativity produce the 
flattening in the correction between the QPO frequency i^qpo = t^K{fm) and the mass accretion 
rate M. But z/qpq becomes truly independent of M only when rin approaches r^so or the stellar 
radius R. It has been suggested (see §1) that the M-independent QPO frequency corresponds to 
the Kepler frequency at riso- But it is also possible that the inner disk radius reaches the stellar 
surface, which is outside the ISO, as M increases. We note that observationally it is difficult to 
distinguish the flattening of i^qpo and a true plateau. It is not clear that the flattening feature 
at i^QPO ~ 1100 Hz observed in 4U 1820-30 (Zhang et al. 1998) corresponds the maximum QPO 
frequency, but we shall assume it does and explore the consequences. 

The maximum QPO frequency, I'ma.x, is given by the orbital frequency at the larger of riso and 
R. To linear order in z^^ (the spin frequency), the ISO is located at riso = (6GM/c^)(l — 0.544a), 
and the orbital frequency at ISO is 

1571 

lyRinso) = -r^(l + 0.748a) Hz, (21) 
with the dimensionless spin parameter 

a^0m9^(^^), (22) 
Mi.4 V300Hzy ' ^ ^ 

where we have adopted / = {2/5)k,AIR'^ for the moment of inertia of the neutron star, with 
K ~ 0.815 (appropriate for a n = 0.5 polytrope). The orbital frequency at the stellar surface can 
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be written, to linear order in z^^, as 



uk{R) = 2169 M^^^ R^o 1 - 0-094 



a 
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Mia\ 



Rio J. 



Hz. 



(23) 



Note that in the above equations, R refers to the equatorial radius of the (spinning) neutron star, 
which is related to the radius, Rq, of the corresponding nonrotating star by: 



where we have again adopted the numerical parameters appropriate for a n = 0.5 polytrope (Lai 
et al. 1994). One may appeal to numerical calculations (e.g.. Miller, Lamb Sz Cook 1998) for more 
accurate results, but the approximate expressions given above are adequate. 

Figure 2 shows the contours of constant i^max = min[z^_ft:(nso)5 i^k{R)] in the M-Rq plane. For 
large M and small Rq, the contours are specified by i^k (fiso), while for larger Rq and small M, 
the contours are specified by vk{R)- We see that to obtain the maximum QPO frequency of order 
1100 — 1200 Hz, one can either have a M ^ 2Mq neutron star (with i?o ^ 16 Km), or have a 
M ~ 1.4M0 neutron star with i?o — 14 — 15 km. Here we focus on the latter interpretation, in 
which the accretion disk terminates at the stellar surface before reaching the ISO. A boundary 
layer forms in which the angular velocity of the accreting gas changes from near the Keplerian 
value (at the outer edge of the boundary layer) to the stellar rotation rate. Depending on the 
thickness of the boundary layer, the inferred the NS radius may be somewhat smaller. Moreover, 
the peak rotation frequency may be below vk{R), which would also allow smaller values of Rq. 

In addition to avoiding a large neutron star mass (see §2), the identification of fmax with 
the Kepler frequency near the stellar surface may allow a plausible explanation of the observed 
correlation between the QPO amplitude and the X-ray flux. While the mechanism of producing 
X-ray modulation in a kHz QPO is uncertain, in many models (e.g.. Miller et al. 1998; see also 
Kluzniak et al. 1990) the existence of a supersonic "accretion gap" between the stellar surface and 
the accretion disk is crucial for generating the observed the X-ray modulation. If we interperate 
t'max as the Kepler frequency at the ISO, which is always outside the stellar surface, then the 
"accretion gap" always exists, and there is no qualitative change in the flow behavior as the inner 
disk approaches ISO. It is therefore difficult to explain why the QPO amplitude decreases and 
eventually vanishes as the X-ray flux increases. The situation is different if t'max = ^k{R), since 
the gap disappears when the mass accretion rate becomes sufficiently large. At small M there is 
a gap (induced by a combination of magnetic and GR effects) between the inner edge of the disk 
and the stellar surface. Since the impact velocity of the gas blob at the stellar surface is larger for 
a wider accretion gap, we expect the modulation amplitude to be larger for small accretion rates 
0. As M increases, the inner disk edge approaches the stellar surface, and we expect the QPO 
amplitude to decrease. The maximum QPO frequency signifies the closing of the accretion gap 



^When M is too low (for a given Bo) so that nn is far away from the steUar surface, the accreting gas can be 




(24) 
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and the formation of a boundary layer. Since there is no supersonic flow in this case, one might 
expect the QPO ampHtude to vanish. In addition, there may be changes in the spectral properties 
of the system as the gap closes. 

The large neutron star radius (15 km for a IAMq star) required if t'max = '^k{R) is only 
allowed for a handful of very stiff nuclear equations of state (sec Fig. 2); most recent microscopic 
calculations give Rq ~ 10 km (e.g., Wiringa et al. 1988). Is such a large radius consistent with 
observations? No neutron star radii are known with the accuracy that has been achieved for 
numerous neutron star mass determinations, but several methods have been tried: 

1. Observations of X-ray bursts have been used to determine empirical M — R relations, 
but these are hampered by the need for model-dependent assumptions regarding the total 
luminosity and its time history, anisotropy of the emission, radiated spectrum and surface 
composition, even when the source distance is known (e.g. van Paradijs et al. 1990, Lewin, 
van Paradijs &i Taam 1995). 

2. X-ray and optical observations of the (apparently nonrotating) isolated neutron star RX 
J185635-3754 (Walter, Wolk & Neuhauser 1996, Walter & Matthews 1997), combined with 
limits on the source distance, D, imply a blackbody radius R{1 + z) < 14(D/130pc) km, 
where z is the surface redshift of the star. 

3. Ray tracing and lightbending may be used to derive limits on R/M for periodically 
modulated X-ray emission. For two isolated neutron stars (PSR B1929-)-10 and B0950-I-08; 
Yancopoulos, Hamilton & Helfand 1994, Wang & Halpern 1997) and one millisecond pulsar 
(J0437-4715; Zavlin and Pavlov 1997, Pavlov & Zavlin 1998), the results arc broadly 
consistent with Re? /2GM ~ 2.0 — 2.5, but the results depend on geometry (angles between 
rotation and magnetic axes, and rotation axis and the line of sight) as well as on the spectrum 
and (energy-dependent) anisotropy of the polar cap emission. The rather large observed 
pulsed fractions appear to rule out two polar cap hot spots unless Rx? /2GM is rather large 
(e.g. ~ 4.3 for PSR B1929+10; Wang & Halpern 1997). Similar considerations may prove 
fruitful for periodically modulated flux from X-ray bursts (e.g. Miller & Lamb 1998); the 
pulse fractions observed so far are large, suggesting non-compact sources (e.g. Strohmayer 
et al. 1999, who find Rx? /2GM ~ 5 for 4U1636-54, corresponding to an implausibly large 
radius of 21 km for M = I.4M0). 

4. Burderi & King (1998) have argued that requiring the Alfven radius to be intermediate 
between R and the corotation radius, Rco = (GM/ri^)^/^, for the 2.5 ms pulsating source 
SAX J1808.4-3658 (discovered by Wijnands &; van der Klis 1998) implies an upper bound 



channeled out of the disk plane by the magnetic field toward the magnetic poles. The detail of the channeling process 
depends on the magnetic field geometry in the disk (such as the radial pitch angle of the field line) . This may quench 
the kHz QPOs and give rise to X-ray pulsation (as in X-ray pulsars). The pulsating X-ray transient system SAX 
J1808.4-3658 may just be such an example. 
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of R < 13.8(M/M0)V3 km, since the pulsations are detected at the same frequency for 
X-ray count rates spanning an order of magnitude. However, their bound depends on the 
model-dependent assumptions that the count rate is strictly proportional to M and the field 
strength in the disk is dipolar [B oc r~^). (See also Psaltis & Chakrabarty 1999.) 

Taken together, the evidence neither supports nor excludes the possibility that i? ~ 15 km for 
M ~ 1.4Mq (or Rc^ /2GM ~ 3.6) definitively, although most of the estimates listed above favor 
more compact models {R ~ 10 km for M ~ 1.4Mq) nominally. 



4.2. QPOs from r > Hn? 

QPOs are identified in the Fourier spectra of photon counts from X-ray sources, so it may 
be that most of the spectral power comes from radii outside nn, possibly from the disk radius at 
which the differential photon emission rate is maximum. For example, if the QPO arises from 
a radius r = (1 + A)rin, then z/qpq = (1 + A)~^/^zvA'(nn)- As 6GM/c^, the ISO in the 

slow-rotation limit, z/qpq 2200 Hz/(M/M0)(l + A)'^/^, so observations that give fmax — 1060Hz 
asymptotically may actually require M(l + A)3/2 ^ 2.IM0, or 1 + A 1.3 if M IAMq. rather 
than M ~ 2.1M0. 

To obtain a simple realization of this idea, consider a Shakura-Sunyaev (1973) disk, for which 
the emitted flux from one face is 

4 SGMMfjr) 
F{r) = asBT, (r) = ^—3 ; (25) 

in the Newtonian limit (which we shall employ here for giving a simplified illustration). The 
function /(r) = 1 — P\/ri^/r, where f3 < 1 parametrizes the rate of accretion of angular momentum 
from the disk onto the star relative to M\/ GMr\^ (e.g. Shapiro & Teukolsky 1983. eq. [14.5.17]; 
see also Frank, King &: Raine 1992, §5.3); if "imperfect" fluid stresses vanish at rin, then /3 = 1 
(as in black hole accretion; see Page & Thorne 1974, Novikov & Thorne 1973). If the color 
temperature of the emission equals the effective temperature Te(r), then the "bolometric flux" 
of photons is ~ F{r) /kT(,{r) at radius r, and the rate at which photons are emitted from radii 
between r and r + dr \s of order 

27rrF(r) M^/^[/(r)]3/4 

^TUO ■ ^ ^ 



Differentiating equation (Pq) implies a maximum emission rate at \/r/r\^ = 1.3/3, consistent with 
f > Tin provided that P > 0.77. Assuming that the QPO frequency is the Kepler frequency at the 
radius of peak (bolometric) photon emission, 

_ l^Kinn) 1000 Hz 
^QPO - (LW " /33M/M0 ' ^^^^ 
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where the limiting result is for rin — > QGM/c^ . In order for the maximum value of z^qpo to be 
z^max 1060 Hz, we require M = 1.4Mo/(/3/0.88)3. 

Real disk emission profiles for small nn, and the determination of fqpo, are not this simple 
for several reasons. A detailed calculation of the X-ray spectrum is needed, since the QPOs are 
found for counts in particular energy bands; the bolomctric count rate is not a good approximation 
in general. (But note that Comptonization by hot coronal gas above the disk conserves photon 
number, so the approximation may be better than it appears at first sight.) In particular, the 
color temperature is not usually the same as the effective temperature, since electron scattering 
is the dominant opacity at relevant disk radii. The composition of the disk is also important; at 
low enough M, the disk will be matter-dominated, but at larger M, radiation-dominated. (Less 
important, but still significant, is the dependence of opacity on the element abundances in the 
accreting gas.) In addition, relativistic effects alter /(r) (e.g. Page & Thorne 1974, Novikov & 
Thorne 1973), and hence fQpo. Moreover, the angular momentum carried away by photons may 
not be insignificant once r-m approaches the ISO (Page &; Thorne 1974, Epstein 1985). Instabilities 
associated with the transition from matter to radiation domination (Lightman & Eardley 1974) 
or the inner boundary layer (e.g. Epstein 1985) might also play a role in determining i^qpo- 
These and other issues associated with the termination of disks at nn and QPOs will be explored 
more fully elsewhere. However, the simplified example presented here indicates that i^qpo might 
plausibly arise from r > rin. 

5. Conclusion 

In this paper we have presented a phenomenological model of the inner region of the accretion 
disk for weakly magnetized neutron stars such as those in LMXBs. A notable feature of these 
systems is that both magnetic field and general relativity are important in determining the inner 
disk radius. Our result suggests that the combined effects of a steepening magnetic field - which 
is likely for disk accretion onto a neutron star - and general relativity can produce the fiattening 
of the QPO frequency fqpo as the mass accretion rate M increases. If the field steepens fast 
enough with decreasing inner disk radius, fqpo may vary little over a fairly substantial range 
of M at values considerably below the Kepler frequency at the ISO due to general relativity. 
Obscrvationally, the correlation between i^qpo and the RXTE photon count rate has been 
well-established, but the scaling between fQpo and M is ambiguous (Mendez et al. 1998c). 
An observational or phenomenological determination of this scaling would be quite useful in 
constraining the magnetic field structure in LMXBs. 

Currently it is not clear whether the plateau behavior in the QPO frequency has been 
observed. But even if i^qpo ~ 1100— 1200 Hz represents the maximum possible QPO frequency, we 
argue that a massive neutron star (M ^ 2Mq) is not necessaily implied. Instead, a M ~ 1.4M0, 
i2o ~ 14 — 15 km neutron star may be a better solution, and is within the range allowed by some 
nuclear equations of state. If this is the case, the maximum QPO frequency signifies the closing of 
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the accretion gap and the formation of a boundary layer. Alternatively, the QPO frequency might 
be associated with the Kepler frequency at a radius somewhat larger than the inner radius of the 
disk, thus allowing lower mass for the accreting neutron star. In either case, better theoretical and 
phenomenological understanding of the termination of magnetized accretion disks is needed before 
observations of maximal kHz QPOs can be interpreted as purely general relativistic in origin, and 
used to deduce neutron star masses. 
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Fig. 1. — The orbital frequency of at the inner radius of the the disk as a function of the "effective" 
mass accretion rate. The inner disk radius, nn, is obtained by solving eq. ( p^ ) (light solid line) or 
([l5|) (dashed solid line), corresponding to different magnetic field structure. The heavy solid lines 
incorporate the effect of general relativity based on eq. (|l^) (heavy solid line) or (|20|) (heavy dashed 
line). 
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Fig. 2. — Constraints on the mass-radius (M — Rq) relation of neutron star from the maximum 
orbital frequency I'max outside the star. Each closed curve shows the fmax =constant contour in the 
Af — Rq plane, with the upper boundary i^max = 1100 Hz, and the lower boundary fmax = 1200 Hz. 
(the solid heavy lines correspond to spin frequency = 300 Hz, and the heavy dashed lines i^g = 0). 
Note that in the case of rotating neutron star, Rq is the radius of corresponding nonrotating stellar 
model with the same mass. The light solid curves depict two representative equations of state: TI 
is the very stiff tensor interaction model of Pandharipande Sz Smith (1975), and UVM+UVH is 
from model of Wiringa et al. (1988). 



